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Collisionless spin dynamics in a magnetic field gradient 
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We study the collisionless spin dynamics of a harmonically trapped Fermi gas in a magnetic held 
gradient. In the absence of interactions, the system evolution is periodic: the magnetization develops 
twists, which evolve into a longitndinal polarization. Recurrences follow. For weak interaction, the 
exchange interactions lead to beats in these oscillations. We present an array of analytic and 
nnmerical techniques for studying this physics. 
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I. INTRODUCTION 

The spins in Fermi gases display a rich range of be¬ 
haviors: Both collisionless and collision dominated spin 
waves have been studied in a number of recent experi¬ 
ments m- Here we present a simple model of the col¬ 
lisionless regime, and study the dynamics. 

Collisionless spin waves were first explored in ^He film, 
then in spin polarized hydrogen [SHlo]. However, some 
of these early works are hard to observe in these atomic 
systems. Unlike these early studies, cold atom experi¬ 
ments allow one to directly image the spin waves mm- 
Such images complement the spectroscopic and transport 
probes used to understand ^He and hydrogen. 

While our focus is the weakly interacting regime, much 
recent work has explored the collision dominated strongly 
interacting limit. One of the problems is to understand 
the spin diffusive behavior in this regime. For example, it 
is conjectured that quantum mechanics forbids the spin 
diffusion constant from exceeding ^ h/m, where h is the 
reduced Planck’s constant and m is the atomic mass. 
This bound appears to be consistent with the recent the¬ 
oretical and experiment results [THSl [TSUTT] . 

Our model is largely inspired by the experiments of 
Bardon et al. in Toronto [S]. There a two-component 
Fermi gas is placed in a cigar shape trap. The gas is pre¬ 
pared with a uniform magnetization in the x direction. A 
magnetic field is applied in the z direction. The strength 
of this field varies linearly with z. The experimentalists 
let the system evolve for some time then interrogate it, 
typically with a spin-echo protocol [miis]. For simplic¬ 
ity, we do not model the spin-echo, but simply study how 
the magnetization evolves. Since all dynamics are one¬ 
dimensional (ID), and the Hamiltonian can be integrated 
out in the other two dimensions, we will focus on the ID 
model instead. 

This paper is organized as follows. First we consider 
the non-interacting case, where we can analytically de- 
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scribe the dynamics. We show that the transverse mag¬ 
netization oscillates, and explain this behavior in terms 
of transverse and longitudinal spin rotations. Then we 
study weak interaction. We find beats which can be 
attributed to the exchange interaction. In the end, we 
discuss the transition from collisionless to the collisional 
limit and give our summary. 

II. NON-INTERACTING CASE 

We consider ID pseudospin-1/2 Fermi gases in a mag¬ 
netic field gradient along the x axis. In the absence of 
interactions, we can consider the dimensionless single¬ 
particle Hamiltonian 

Ho = -f + -+aXx, (1) 

where tr = ± for the up and down spin states, with re¬ 
spect to the z axis. Here x = xj^hjivrvjS), where x is 
the position of the particle and ^/hj (moj) is the char¬ 
acteristic length of the harmonic oscillator, with uj the 
trap frequency. The magnetic field gradient is repre¬ 
sented by dimensionless A. The actual magnetic gradient 
is A = Xujy/hmjj/pL, where pL is the magnetic moment. 
The dimensionless time t = Lot, where t is the real time. 
Hereafter we will use these dimensionless quantities. The 
many-body state will consist of a Slater determinant of 
single-particle states tpna-iXjt), each of which evolve via 
the time-dependent Schrodinger equation, idttpn = Hif^- 
The index n runs from 1 to N, where N is the number 
of particles. 

We envision preparing the system so that at time t = 0 
the n’th state of the Slater determinant has spin com¬ 
ponents V’nt = 4>nix)lV^ and ipni = 4'nix)lV2, where 
4>n{x) is the n’th eigenstate of the simple harmonic os¬ 
cillator. Such an initial condition is prepared by first 
polarizing the system in the f state, and cooling it to the 
ground state. One then applies a 7r/2 pulse, rotating the 
spin from the z direction to the x direction. 

To find the subsequent dynamics, we make the follow¬ 
ing ansatz: 

V’n-T {x,t) = (/„ (x - (TXo (t)) . (2) 
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This ansatz describes the exact dynamics in the ab¬ 
sence of interactions. Substituting this ansatz into the 
Schrodinger equation yields a set of ordinary differential 
equations for xq, v and 

Xo(t) + xo(t) + A = 0, (3) 

xg(t)/2-u^(t)/2-<(t)-if„ = 0. (4) 

These equations are readily integrated to yield 

xo (t) = A (cost — 1) ,v (t) =—Asint, (5) 

(t) = - - -A smt- i^n+-- —jt. (6) 

The local magnetization in the x, y plane can be ex¬ 
pressed as a complex number m(x,t) = mx + imy. 
This complex magnetization can be written as m(x,t) = 

J2n'^ri(x,t) with 

mn(x,t) = i^^^(x,t)il;nt(x,t) 

= (t)n{x- Xo (t)) (j>n (x + Xq (t)) ( 7 ) 

We give the physical picture leading to magnetization 
dynamics in this non-interacting case in Fig. As illus¬ 
trated in Fig. [^a), the center of the up-spin and down- 
spin clouds become separated by a distance 2xo(t). This 
reduces the overlap between the clouds, and the magni¬ 
tude of the local transverse magnetization. As is evident 
in Eq. Q, the clouds also develop a relative phase profile 
exp(2z;^. As shown in Fig.[^b), this phase factor can be 
interpreted as a spin precession term. While it does not 
change the magnitude of the local polarization, it makes 
the direction depend on position. Thus it does reduce the 
total transverse magnetization M(t) = f m(x,t) dx. By 
symmetry M(t) is always real, implying the net trans¬ 
verse polarization is always in the x direction. By con¬ 
struction M(t) < N is dimensionless. 


from Eq. Q, we see clear undamped oscillations. We 
also show the contribution from the n = 0 and n = 1 
states. The lower energy states dominate the dynam¬ 
ics, as the maximum displacement \xo(t = 7r)| = 2A is a 
greater fraction of their width. 


M{t) 



FIG. 2: Time evolution of average transverse magnetization 
per particle M(t)/N (solid) for A = 21 particles. Here t = cut 
is the reduced time, where t is the real time and lu is the trap 
frequency. Also shown is the average magnetization of the 
particles for the two lowest states: n = 0 (dotted) and n = 1 
(dashed). 


III. WEAK INTERACTION 

In this section, we consider the influence of weak in¬ 
teraction on the system. We make a time-dependent 
Hartree-Fock ansatz, assuming that at all times the sys¬ 
tem is described by a Slater determinant of N single¬ 
particle states. The equations of motion can be derived 
by minimizing the action 



FIG. 1: Illustration of the processes which lead to magneti¬ 
zation dynamics of a Fermi gas in a magnetic field gradient, 
(a) Longitudinal spin dynamics: The up-spin and down-spin 
atoms move relative to each other, changing their overlap, and 
hence the transverse magnetization, (b) Transverse spin dy¬ 
namics: The spins precess in the transverse plane, at speeds 
which depend on the location of the atoms. This leads to an 
inhomogeneous texture, whose average magnetization is re¬ 
duced. The single arrow labels the spin direction, while the 
double one represents the movement of the spin. 


S = J ^*(idt-H)'i>dtd^x, (8) 

where '^'(xi, X 2 , ■'' ^XN,t) represents the many-body 
wave function. The action can be decomposed as a sum 
of a non-interacting and an interacting part S = S'o + S'/. 
Within the Hartree-Fock ansatz, 

•5o = + 

m,n 

with g = g/(hojy/H/(muj)) represents the reduced inter¬ 
action strength with g the physical one. Minimizing the 
action with respect to the wave functions '0* ^ (x) yields 
equations of motion 


— aAxj tpn,( 7 dxdt, 


— “b ^ ^ (|V^m,cr| V^n,cr 

m 

''Pm,a^'rn,s'4^n,s') ■ 


Figure[^shows the evolution of the transverse magneti¬ 
zation for TV = 21 particles, taking A = 0.1. As expected 


idttpn,(T 


(9) 
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As argued in [^, in the limit of slow spatial and tem¬ 
poral dynamics, these coupled equations are equivalent 
to the collisionless Boltzmann equation. For moderate 
N, we can numerically integrate these equations. The 
exchange interaction acts like a spatially inhomogeneous 
transverse magnetic field, and scrambles the spins. The 
simple oscillations seen in Fig. [^develop more structure, 
exhibiting quantum beats as seen in Fig. 

These beats can be qualitatively understood by consid¬ 
ering a single wave function '4’n,a{x) in the Slater deter¬ 
minant, and treating the other wavefunctions as a static 
homogeneous background. Within this approximation, 
Eq. § becomes a similar Jaynes-Cummings model. 


idt 


V 



where 


- d‘^ 

H = —^ ^ z -f ge<Jx, 


( 10 ) 


( 11 ) 


with the effective exchange interaction strength. The 
subsequent mathematics is simpler if we rotate spin space 
by 90°. This is done by introducing a transformation 
matrix 

and defining 



H' 



RHR. 


Equation (10) then becomes 


idt 


'>Pn- 


= H' 


V'n.- 

V'n.-f 


(13) 

(14) 


(15) 


where 


/ {n+\)+ge +a) \ 

V {n+\)-g^ j 


(16) 


with a = [x + dx)/ -s/S, = {x — dx)/ V^- In this rotated 

vector space our initial state is = 0) = 

= 0 ) = 0 . 

This Hamiltonian can be understood from the level di¬ 
agram in Fig. For zero interaction and magnetic field 
gradient, the system is a simply harmonic oscillator and 
has energy levels at En = n + 1/2, which is the dashed 
lines in Fig. The exchange coupling ge shifts the eigen¬ 
state energy of —>■ and ^ spins, shown as the solid lines 
and energy En ± ge- The magnetic field gradient flips 
the spins and at the same time changes the vibrational 


level by one. For small A, one can truncate the time dy¬ 
namics to a three level system: the initial state, labeled 
by a circle, and the two states which are coupled to it. 



FIG. 3: Level structure of the simplified model in Eq. (151, 
used to describe the time evolution of a single wavefunction 
'>pn,a{x) in the Slater determinant. Dashed lines show the 
eigen-energies when the field gradient A and effective exchange 
field Qe are both zero. The exchange field splits the degeneracy 
between —>■ and states, yielding the solid lines. En = n -\- 
1/2. The gradient A couples levels as shown by the arrows. At 
time t = 0, 'tpn,^{x) = 4in{x) and = 0, corresponding 

to the state marked with a circle. 


This yields two oscillation frequencies i/± = lit 2ge, and 
to second order in A the contribution to the magneti¬ 
zation from the n’th term in the Slater determinant is 


m 


= 1 cos(l - 2ge)t + 


(l_2gD^ CUO^J- - z,ye;i- w (i+ 2 gJ 2 COs(l -|- 2ge)t. 
The beats in Fig. are qualitatively consistent with this 
scenario, but this simplified model is not able to quanti¬ 
tatively describe the evolution of the magnetization. A 
more sophisticated, and quantitatively accurate, ansatz 
involves taking the n’th wavefunction in the Slater deter¬ 
minant as having the form 


n ,(7 , t'j A^nit) 

-I- crH„(t)(/i„+i(a:) -I- aCn{t)4>n-i{x). (17) 


Minimizing the action with this ansatz yields equations 
of motion for An, Bn, and Cn, 


/ idt — En —XXn —XXn-1 
I XXn idt En-\-l 
y XXn—1 idt En—1 


1 

^ ( Anin I Bn 

™ J \ Cn 


= 0 , 


where A„ = f dx(j)n{x) 4 >n+i(,x). The nonlinear term is 
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A„™= + 

\ + Cmln-1 ) 


m) 



FIG. 4: (Color online) Time evolution of average trans¬ 
verse magnetization for reduced interaction strength g = 
g/{hujy/h/{mLo)) = 0.1. The particle number and magnetic 
gradient is the same as the non-interacting case and t = ujt 
is the reduced time with t the real time and oj the trap fre¬ 
quency. Solid black line: Simplified variational approxima¬ 
tion. Dashed orange line: Full Hartree-Fock calculation. 

where 

(3™ = y dx4>l^{x)(l)n+l(,x)(j)n-l{x), 

7 ™ = y dx4>m(.x)(l)m+l{x)4>n{x)4>n+l{x)- 

The above equations can be solved iteratively and we get 
the net magnetization 

M{t) = ^ (|A„(t)|2 - |i?„(t)|2 _ |C„(t)p) . (19) 

n 

We plot the magnetization for system N = 21 and 
A = 0.1 in Fig. The reduced interaction is chosen 
as g = 0.1. We can see the weak interaction develops 
beats in the magnetization, which is expected from our 
previous approximation. These differential equations are 
much easier to integrate than the Hartree-Fock equations, 
and as seen in Fig. this simple ansatz captures most 
of the relevant physics. 

Finally, we give the condition when the actual three- 
dimensional (3D) experiment goes from collisionless to 
collisional limit. As argued in this transition is 

determined by F/w, where F = nav is the two-body 
elastic scattering rate per particle with n the mean den¬ 
sity of the system, a the scattering cross section, and 
V the relative velocity. If the scattering event seldom 


I__ M 

happens at several oscillation periods, i.e., F <C w, the 
dynamics is mainly dominated by collisionless processes. 
The system shows spin wave dynamics in this regime. 
Increasing the scattering rate F, it continuously goes 
from collisionless limit to collisional limit. In our sys¬ 
tem n « {2mEp )^^‘^/« 43.8(mw/Ii)^/^/7r^ for 
Ep = 20.5Huj, and v ~ X^/huJ/m. Our ID reduced 
interaction strength is conn ected to the 3D scattering 
length as g = ghjjyJh/irnuj) « 2fi^a^p,/{ral‘^) [^. Here 
I = y/h/(mu}±) is the transverse width and for a isotropic 
harmonic potential we have a^pi = gy/K/JrruJ) j2. Thus 
we have F = Airna^jjV ~ 43.8(7^Aw/tt. So for a 
strongl y interacti ng collisional limit, one has to fulfill 
g Y^7r/(43.8A) « 0.8 for reduced magnetic gradient 
A = 0.1. Beyond this point, the system is dominated 
by collisional processes and the physics of spin diffusion 
begin to take part. 

IV. SUMMARY 

In summary, we give an analysis of the collisionless 
spin dynamics of a weakly interacting Fermi gas in a 
magnetic field gradient. In the absence of interactions 
the net transverse magnetization oscillates with the trap 
frequency. Weak interactions lead to beats. We present 
simple models which explain these dynamics. 

For a related experiment, the magnetization dynam¬ 
ics can be directly observed using the techniques in Ref. 
[5], and the beats behavior clearly tells the low energy 
spin wave dynamics of the system. Here even though 
our calculation is purely ID, we believe the beats behav¬ 
ior is also preserved in the 3D system. This can also 
be seen from our mean-field argument, where the other 
dimensions will only change the eigenenergy of the har¬ 
monic oscillator. However the other two dimensions do 
contribute to the system. They will influence the local 
interaction strength and then shift the beats. Also in 
the experiment one needs to take into account the effect 
of temperature. For a non-zero temperature, the atoms 
will spread to higher single-particle states and the spin 
dynamics will be enhanced. 
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